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Abstract 

We analyse within the "flavoured" leptogenesis scenario of baryon asymmetry 
generation, the interplay of the "low energy" CP-violation, originating from the 
PMNS neutrino mixing matrix U, and the "high energy" CP-violation which can 
be present in the matrix of neutrino Yukawa couplings, A, and can manifest itself 
only in "high" energy scale processes. The type I see-saw model with three heavy 
right-handed Majorana neutrinos having hierarchical spectrum is considered. The 
"orthogonal" parametrisation of the matrix of neutrino Yukawa couplings, which 
involves a complex orthogonal matrix R, is employed. In this approach the matrix 
R is the source of "high energy" CP-violation. Results for normal hierarchical (NH) 
and inverted hierarchical (IH) light neutrino mass spectrum are derived in the case 
of decoupling of the heaviest RH Majorana neutrino. It is shown that taking into 
account the contribution to Yg due to the CP-violating phases in the neutrino 
mixing matrix U can change drastically the predictions for Yg, obtained assuming 
only "high energy" CP-violation from the iZ-matrix is operative in leptogenesis. 
In the case of IH spectrum, in particular, there exist significant regions in the 
corresponding parameter space where the purely "high energy" contribution in Yb 
plays a subdominant role in the production of baryon asymmetry compatible with 
the observations. 
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1 Introduction 



In the present article we investigate further the possible connection between leptogene- 
sis [1, 2] (see also, e.g. [3, 4]) and the low energy CP- violation in the lepton (neutrino) 
sector (for earlier discussions see, e.g. [5, 6, 7, 8] and the references quoted therein). It 
was shown recently in [9] that the CP-violation necessary for the generation of the ob- 
served baryon asymmetry of the Universe in the thermal leptogenesis scenario can be 
due exclusively to the Dirac and/or Majorana CP- violating phases in the Pontecorvo- 
Maki-Nakagawa-Sakata (PMNS) neutrino mixing matrix [11], and thus can be directly 
related to the low energy CP-violation in the lepton sector, e.g. in neutrino oscillations, 
etc. (see also [10]). The baryon asymmetry is produced in the regime when the lepton 
flavour effects in leptogenesis [12, 13, 14, 15] are significant ("flavoured" leptogenesis). 
As was realised in [14, 15], the lepton flavour effects can play very important role in the 
leptogenesis mechanism of baryon asymmetry generation. 

In ref. [16], which was stimulated by some of the results obtained in [9], the dependence 
of the baryon asymmetry produced in "flavoured" leptogenesis on the lightest neutrino 
mass, min(mj), j = 1,2,3, when the requisite CP-violation is provided entirely by the 
PMNS matrix f/pMNS? was studied. In certain rather general and physically interesting 
cases the generated baryon asymmetry was found to depend strongly on the value of 
min(mj). For specific values of min(mj), in particular, the asymmetry can be strongly 
enhanced (by a factor of ~ 100 or more) with respect to that predicted in the case 
of min(mj) = 0. This enhancement can make the predicted Yb compatible with the 
observations even when this is not the case for min(mj) = 0. Some aspects of the matter- 
antimatter asymmetry generation in the "flavoured" leptogenesis scenario in the case 
when the relevant CP-violation is due to the Majorana or Dirac CP-violating phases in 
t^PMNS, were investigated also in [17, 18]. 

Here we analyse in detail the interplay in "flavoured" leptogenesis of the "low energy" 
CP-violation, originating from the PMNS neutrino mixing matrix, and the "high energy" 
CP- violation which can be present in the matrix of neutrino Yukawa couplings. A, and 
can manifest itself only at some "high" energy scale. Both types of CP-violation can, in 
general, have their origin in the neutrino Yukawa couplings. The latter, as is well-known, 
are one of the basic ingredients of the see-saw mechanism of neutrino mass generation [19] , 
on which the leptogenesis theory is based. A widely recognised appealing features of 
the see-saw model include i) a natural explanation of the smallness of neutrino masses 
(see, e.g. [20, 21, 22]), and that ii) it allows to relate, through the leptogenesis theory, 
the generation and the smallness of neutrino masses with the generation of the baryon 
(matter-antimatter) asymmetry of the Universe, Yb- 

The minimal scheme in which leptogenesis can be realised is the non-supersymmetric 
version of the type I see-saw model with three (or two) heavy right-handed (RH) Majorana 
neutrinos, Nj, having masses Mj. The matrix of neutrino Yukawa couplings. A, together 
with the Majorana mass matrix of the RH neutrinos, M^, and the matrix of charged 
lepton Yukawa couplings, A'*^^, plays a crucial role both in the see-saw mechanism and in 



1 



leptogenesis. In the basis in which Mr and A'^^ are diagonal ^, the matrix of neutrino 
Yukawa couphngs A is the only source of CP-violation in the lepton sector. Among the 
several possible parametrisations of A (see, e.g. [7]), the orthogonal one, involving a 
complex orthogonal matrix R [23], allows to relate in a rather direct manner the matrix 
A with the neutrino mixing matrix t/pMNS = U: X = (1/w)a/M R^/rriU\ where M and 
m are diagonal matrices formed by the masses Mj > and > of Nj and of the 
light Majorana neutrinos u/^, j, k = 1,2, 3, and v is the vacuum expectation value of the 
Higgs doublet field. This parametrisation proved rather convenient in the analysis [9] 
of the possibility that the CP-violation necessary for a successful leptogenesis could be 
provided by the Majorana and/or Dirac physical phases in the neutrino mixing matrix 
U. It permits to investigate also the combined effect of the CP-violation due to the 
neutrino mixing matrix U and the CP-violation due to the matrix R in the generation of 
the baryon asymmetry in "flavoured" leptogenesis. We will use the terms "low energy" 
and "high energy" for the CP-violation originating respectively from the matrices U and 
R. The PMNS matrix U is present in the weak charged lepton current and can be a 
source of CP-violation in, e.g. neutrino oscillations at "low" energies ^ ~ Mz (see, e.g. 
[26, 27, 28]). The matrix R, as is well-known, does not affect the "low" energy neutrino 
mixing phenomenology. The two matrices U and R are, in general, independent. It should 
be noted, however, that in certain specific cases (of, e.g. symmetries and/or texture zeros) 
of the matrix A of neutrino Yukawa couplings, there can exist a relation between (some 
of) the CP-violating phases in U and (some of) the CP-violating parameters in R (see, 
e.g. [8, 29]). 

The source of the requisite CP-violation in "flavoured" leptogenesis can, in principle, 
be the matrix R, the PMNS matrix U, or both the matrices R and U. The division 
between CP-violation due to the PMNS matrix U and that due to the matrix R at "high" 
energies, e.g. in leptogenesis, needs some clarification. If the matrix R satisfies the 
general CP-invariance constraints (having real or purely imaginary elements [9]), while 
the PMNS matrix, and correspondingly, the matrix of neutrino Yukawa couplings A do 
not satisfy these constraints, we will consider the CP-violation as originating from the 
neutrino mixing matrix U, i.e. from the Dirac and/or Majorana phases in U. If, however, 
the Dirac and Majorana phases in U take CP conserving values, while the matrix R, and 
the Yukawa couplings A do not satisfy the constraints following from the requirement 
of CP-invariance, the CP-violation will manifest itself only in "high" energy phenomena 
(like, e.g. leptogenesis) and will be due to the matrix R. In this case there will be no 
(observable) effects of violation of CP symmetry at "low" energies in phenomena caused 
by the neutrino mixing (neutrino oscillations, neutrinoless double beta {{PP)ou-) decay 
[30, 31], etc.). When neither U nor R satisfy the CP-invariance conditions, both U and 
R will be sources of CP-violation effects at "high" energies . 

In the present article we investigate the possible interplay between the "low energy" 
CP-violation due to the Dirac and/or Majorana CP- violating phases in the PMNS matrix 

^This basis can be chosen without loss of gencraUty and we wiU use it in our further analysis. 
^As is well-known, only the Dirac phase in U can be a source of CP-violation in neutrino oscillations; 
the probabilities of oscillations of flavour neutrinos do not depend on the Majorana phases in U [24, 25]. 
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U, and the "high energy" CP-violation originating from the matrix R, in "flavoured" 
leptogenesis We work within the simplest type I see-saw scenario with three heavy RH 
Majorana neutrinos Nj, j = 1,2,3. The latter are assumed to have a hierarchical mass 
spectrum. Mi -C M2,3. In what concerns the light neutrino masses, we consider two types 
of spectrum allowed by the existing data (see, e.g. [20]), namely, the normal hierarchical 
(NH), mi <C m2 < m^, and the inverted hierarchical (IH), -C rrii < m2, and present 
detailed results for these two spectra. 

There is practically no overlap between the results obtained in our work and those 
found in ref. [15]. One of the objectives of the analysis in [15] was to find out whether 
there can be large differences between the baryon asymmetry of the Universe, predicted in 
the flavoured leptogenesis scenario (the two- or three- flavour regime), and the asymmetry 
derived in unfiavoured leptogenesis (the one-flavour regime). The authors of [15] were not 
interested in (and did not try to provide an answer to) the question of the the relative 
magnitude of, and the interplay between, the contributions to the baryon asymmetry due 
to the "low-energy" CP violating phases in the PMNS matrix and that due to the "high 
energy" CP violating phases in the R matrix, which is the main subject of investigation 
in our work. 

Negative results regarding the possible interplay between the "high energy" and "low 
energy" CP violation in "flavoured" leptogenesis with hierarchical heavy (RH) and light 
Majorana neutrinos were reported in [32]. The results reported in [32] were obtained for 
i) min(mj) = and specific texture zero in the i?-matrix, or ii) for a value of the mass of 
the lightest RH neutrino ^ Mi = 10^° GeV. In case ii) the lightest neutrino mass min(mj) 
was allowed to vary within the interval < min(mj) < 10~^ eV. However, in both these 
cases the contribution in Yb due to the "low energy" CP violation is strongly suppressed. 
The contribution under discussion can be relevant for the production of Yb compatible 
with the observations provided [9] Mi ^ 4 x 10^'^ GeV. In case i) it can be relevant if one 
considers values of Mi ^ 4 x 10^° GeV and of min(mj) ^ 5 x 10~^ eV [16]. 

The lepton flavour effects can be significant in leptogenesis in the case of hierarchical 
spectrum of the heavy Majorana neutrinos, provided the mass of the lightest one Mi satis- 
fies the constraint [14, 15] (see also [18]): Mi ^ 10^^ GeV. In this case the predicted value 
of the baryon asymmetry depends explicitly (i.e. directly) on U and on the CP-violating 
phases in U. Using this fact it was shown in [9] that the observed baryon asymmetry Yb 
can be produced even if the only source of CP-violation in leptogenesis is the Majorana 
and/or Dirac phase (s) in the PMNS matrix ^ U. These results were demonstrated to 
hold both for normal hierarchical (NH) and inverted hierarchical (IH) spectrum of masses 
of the light Majorana neutrinos. They were obtained for CP-conserving elements of the 
orthogonal matrix R. The CP-invariance constraints imply [9] that the matrix R could 
conserve the CP-symmetry if its elements are real or purely imaginary. 

■^Private communication by S. Davidson. We thank S. Davidson for clarifications regarding the analysis 
performed in [32]. 

^The same conclusion was shown to be valid also for quasi-degenerate in mass heavy Majorana neu- 
trinos [9]. 
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Our analysis is performed for negligible renormalisation group (RG) running of rrij 
and of the parameters in the PMNS matrix U from Mz to Mi (see, e.g. [33, 34, 35]) ^. 
Throughout the present work we use the standard parametrisation of the PMNS matrix: 



U 



C12C13 

"■512C23 — Ci2S23Si3e 



S12S23 - Ci2C23Si3e' 



S12C13 
C12C23 ~ Si2S23Sise 



-C12S23 - Si2C23Si3e' 



■523C13 

C23C13 / 



diag(l, e*~, e*~ 



where Cij = cos%, Sij = sinOij, 9ij = [0, vr/2], 6 = [0, 27i] is the Dirac CP-violating (CPV) 
phase and 021 and 031 are the two Majorana CPV phases [24, 36], 021.31 = [0,47r]. All 
our numerical results are obtained for the best fit values of the solar and atmospheric 
neutrino oscillation parameters [37, 38, 39], AttIq, sin^ ^12 and Am^, sin^26'23: 

Am| = Aml^ = 8.0 x 10^^ eV^ sin^ Ou = 0.30, (2) 

\Aml\ = |Am^^32)| = 2.5 x 10"^ eV^ sin^ 2^23 = 1- (3) 

In certain cases ^ the predictions for \Yb\ are very sensitive to the variations of sin^ 2^23 
within its 95% C.L. allowed range: 

0.36 < sin^ ^23 ^ 0.64, 95% C.L. (4) 

We also use the upper limit on the CHOOZ mixing angle 6^13 [41, 37, 40]: 

sin^ 9i3 < 0.035 (0.056) , 95% (99.73%) C.L. . (5) 



2 Baryon Asymmetry from "Low" and "High" En- 
ergy CP- Violation 

Following [9, 16], we perform the analysis in the framework of the simplest type I see-saw 
scenario. It includes the Lagrangian of the Standard Model (SM) with the addition of 

^We have checked by exphcit calculations using the equations describing the RG running of the 
neutrino mixing parameters given in http://www.ph.tum.de/ihratz/AnlyticFormulac/ (sec [34]), that 
the running of the neutrino mixing angles, the Majorana and Dirac phases in the neutrino mixing matrix 
and of the neutrino masses, is negligible in both cases of NH and IH light neutrino mass spectrum studied. 
We have verified, in particular, that for the values of sin ^13 of interest for our discussion, the running 
of the Dirac phase in the case of NH spectrum is so small that can be safely neglected. The running of 
the Dirac phase in the case of IH spectrum with negligible lightest neutrino mass ma, as is well-known, 
is negligible. 

^Using the latest data from the KamLAND and SNO experiments in the global neutrino oscillation 
analysis one obtains somewhat different best fit values of Amg, and |Am^| [40]: Auiq = 7.65 x 10~^ eV^, 
[Am^l = 2.4 X 10"'^ eV-^. The results of our analysis do not change if we use these best fit values. 
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three heavy right-handed Majorana neutrinos Nj (j = 1,2,3) with masses < Mi < 
M2 < M3 and Yukawa couphngs A^;, / = e,fi,T. We will work in the basis in which i) 
the Yukawa couplings for the charged leptons are flavour-diagonal, and ii) the Majorana 
mass term of the RH neutrino fields is also diagonal. The heavy Majorana neutrinos are 
assumed to possess a hierarchical mass spectrum. Mi <^ M2 <S M3. 

In what follows we will use the well-known "orthogonal parametrisation" of the matrix 
of neutrino Yukawa couplings [23]: 

X = -VMRy^U^ , (6) 

V 

where R is, in general, a complex orthogonal matrix, R R^ = R^ R = 1, M and m are 
diagonal matrices formed by the masses of Nj and of the light Majorana neutrinos uj, 
M = Diag(Mi,M2,M3), m = Diag(mi, ms, ms), M,- > 0, mj. > 0, and v = 174 GeV is 
the vacuum expectation value of the Higgs doublet field. In contrast to [9, 16], we shall 
assume that the matrix R is not CP-conserving and thus has complex elements. 

In the case of "hierarchical" heavy Majorana neutrinos Nj, the CP- violating asymme- 
tries, relevant for leptogenesis, are generated in out-of-equilibrium decays of the lightest 
one Ni. The asymmetry in the lepton flavour / (lepton charge Li) is given by [14, 15]: 

1/2 3/2, 



3Mi ^^[Ejk^j W UljUikRijRik^ 
- ~T^Z3 TTT^ • (7) 



Obviously, the total asymmetry (eg + + Ct) does not depend on the neutrino mixing 
matrix U (E. Nardi et al., ref. [14]). 

We shall assume that the baryon asymmetry is produced in the "two-flavour" regime 
in leptogenesis [14, 15]. This regime is realised at temperatures 10^ GeV < T ~ Mi < 
10^^ GeV. For T ~ Mi from the indicated interval, the Boltzmann evolution of the 
asymmetry in the r— flavour (lepton charge of the Universe) is distinguishable from 
the evolution of the (e + /i)— flavour asymmetry eg + e^. In the two-flavour regime, the 
baryon asymmetry ^ predicted in the case of interest is given by [15]: 

_ 12 / /417 ^ \ /390 ^ 

Here g^, = 217/2 is the number of relativistic degrees of freedom, €2 = eg -|- e^, m2 = 
rhg + m^, fhi, I = e,fi,T, is the "wash-out mass parameter" for the asymmetry in the 
lepton flavour / [14, 15], 

2 



mi 



l = e,n,T, (8) 



^It is well known that, in thermal leptogenesis with hierarchical heavy Majorana neutrinos, the con- 
tribution of the next to lightest one, N2, to the baryon asymmetry can be significant [42]. We limited 
our analysis to the conventional TVi dominated scenario because we didn't want to assume particular 
constraints on the neutrino Yukawa couplings, such that, for example, can imply no iVi washout effects 
during N2 leptogenesis. 

®The expression we give is of the baryon asymmetry normalised to the entropy density, see, e.g. [9]. 
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and ?7(390mT-/589) = 77(0.66771,-) and ?7(417m2/589) = 77(0. 71m2) are the efficiency factors 
for generation of the asymmetries and 62- The efficiency factors are well approximated 
by the expression [15]: 



8.25 X 10-W / X ^ ^-^^ 



For the observed value of the baryon asymmetry Y^^^, to be reproduced by leptogenesis, 
we will use in our analysis Y^^* = 8.6 x 10~^^ and the "conservative" interval 8.0 x 10~^^ < 
Y^'" < 9.2 X 10-". 

For complex Rij of interest it proves convenient to cast the asymmetries ei, I = e, fi, r, 
in the form: 



= ^ip 12 S ^ I r I sin 2(^1^ + V V ^m/^mp I I 

(10) 

X [{nip - nip) cos((/9/3p) Im (f/;^ f/fp) + {nip + m^) sm{(ppp)Re (U^^Uep)] | , 
where we have used 

Rij = \Rij\e"^^' , ipij = (pii + (pij. (11) 

The ffist term in the curly brackets in eq. (10) represents the contribution to from the 
"high energy" CP-violation, originating entirely from the matrix i?, while the terms in 
the square brackets are "mixed" , i.e. they are due both to the "low" and "high" energy 
CP-violation, generated by the neutrino mixing matrix U and by the matrix R. In what 
follows we will call the CP-violating phases (pij associated with the matrix R, "i?-phases" 
or "high energy" phases. Obviously, if (pij = kj'7i/2, kj = 0, 1,2, j = 1,2,3, the "high 
energy" term is zero, while the "mixed" term reduces to a "low energy" term in the sense 
that, with exception of very special cases (see [9]), the only source of CP- violation in 
leptogenesis will be the PMNS matrix U. The expression (10) for ee implies that, since 
Rij satisfy the orthogonality condition Rf^ + RI2 + -R13 = 1, we can have 7^ only if 
at least two of the three elements Rij of the first row of R are different from zero. We 
will be interested primarily in the effect the "high energy" and the "mixed" terms have 
on the predicted value of Yb- 



3 Normal Hierarchical Neutrino Mass Spectrum 

In the case of NH light neutrino mass spectrum we have mi -C m2 < m^, and consequently, 
7712 = A/Am| = 9 X 10"^ eV, Tng = ^/Anl\ = 5 x 10"^ eV. We shall investigate here the 
case of negligibly small lightest neutrino mass 777.1. More precisely, we will assume that 
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the terms proportional to (mi)"', n = 1/2; 1; 2, in the expressions for e,-, 62, and m2 
are significantly smaller than the terms proportional to (^2,3)"', and can be neglected. In 
what follows we set mi = for simplicity. The asymmetry e,- in this case takes the form: 



3 Ml v^Ami r f Am 



1-^12^ \ Ut2\'^ sin 2^912 



|i?i3|' It/.sP sin 2(^13 + ( ^;;^J 1^121 l^isl 



AmQ 



Am| 
Am\ 



sin(v9i2 + V5i3) Re ([/;2 t^rs) 



cos(v9i2 + V5i3) Im (f/^a f^rs) 

(12) 



It is easy to show, taking into account the unitarity of the matrix U, that the expression 
for the CP-asymmetry €2 can be simply obtained from the expression for e,- : 

e2 = ee+e^ = e,(|f/.fe|' ^ 1 - \Urk\^ , f/;2f/.3 ^ -K^Urs) , = 2, 3 . (13) 

The first term in the brackets in eq. (12) is suppressed by the factor Am^/ Am\ = 
0.03. A more detailed study shows that it always plays a subdominant role in the gener- 
ation of baryon asymmetry compatible with observations and can be safely neglected. 

It should be clear from the expressions (12) and (13) for the asymmetries and €2 that 
the CP-violation due to the PMNS matrix U can play a significant role in leptogenesis 
only if the "mixed" term proportional to I-R12-R13I in eq. (12) is comparable in magnitude, 
or exceeds, the "high energy" term oc | i?i3p| [7^3 P sin 2*^13. The latter will not give a 
contribution to the asymmetries e,- and €2 if sin2(^i3 = 0, i.e. if -R13 is real or purely 
imaginary [9]. 

The elements of the matrix R are constrained by the orthogonality condition: Rf^ + 
R^2 + -^13 = 1- the case of "small" lightest neutrino mass mi under consideration, the 
Rii element does not appear in the expressions for Sr, €2, m-r and m2, which are relevant 
for the calculation of the baryon asymmetry Yb- We will analyse in what follows the 
possibility of relatively small |-Rii|, so that the term R^^ in the orthogonality condition 
can be neglected. This is realised if l-Rup ^ min(l, |i?i2p| sin 2i^i2|). Under this condition 
we can set i?ii = in all our further considerations. Let us note that we get i?ii = 0, 
e.g. in the case of decoupling of the heaviest RH Majorana neutrino [8, 43], leading 
effectively to the so-called "3 x 2" see-saw model ^ [45]. 

For negligible |-Riip, the orthogonality condition for the elements of R can be written 
in terms of two equations involving the absolute values and the phases of i?i2 and -R13: 



^In the context of "flavoured" leptogenesis the case of A^3 decouphng was discussed earher in [15, 9]. 
However, our analysis practically does not overlap with the analyses performed in [15, 9]. Let us note 
also that in the SUSY version of the "3 x 2" see-saw model it can be experimentally feasible, in principle, 
to reconstruct the i?-matrix, as discussed in [44]. If achieved, such a reconstruction could shed light on 
the origin of the CP violation in leptogenesis. 
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(14) 
(15) 

Obviously, we should have sgn(sin 2(^12) = —sgn(sin 2(^13). Using these equations we can 
express the phases (pyi and (pn in terms of I-R12P and |-Ri3p: 





COS 2(^12 + 




COS 2(/)i3 


= 1, 




sin 2(^12 + 


■ -Ri3p 


sin 2(^13 


= 0. 



cos 2(^12 = ^ — ^— ■, sin 2(^12 = ±v^l - cos^ 2v9i2 , (16) 



1 + 


-R12 ^ ^ 






2|i?12P 




1 - 


ii?i2r + 





cos2v9i3 = ^ — ^— , sin2<^i3 = =fV'1 - cos^ 2v?i3 • (17) 

The fact that — 1 < cos2<^i2(i3) < 1 leads to the following conditions: 

(i+ii?i2r)'> i/2i3i', (1 - ii?i2r)' < ii2i3h (18) 

(l + |i?i3r)'> |i?i2|', (1- 1^131')' <|i?12|'. (19) 
Alternatively, one can express I-R12P and I-R13P as functions of the phases: 

2 sin 2(^13 



\R 



12 



\R 



13 



sin2(v9i3 - ifvi) ' 

sin 2(^12 
sin2(<^i3 - (^12) 



(20) 



We will consider values of the phases </?i2 and (y9i3 from the interval [0,27r]. The positivity 
of |i?i2p and I-R13P implies the following constraints on the allowed ranges of (pyi and (^13: 

k-K < <pi3 < (2A; + l)f , (^13 - f - k'n < (^12 < {k - k')n ; (21) 
(2A: + l)f < (^13 < (A: + l)7r, (A; - A;')7r < (^12 < ^13 - f - , (22) 

where A; = 0, 1, 2, 3 and k' = 0, ±1, ±2, ±3. 

We will be interested in the case when the "mixed" term oc I-R12-R13I in the expres- 
sion (12) for the CP-asymmetry is sufficiently large and gives either a dominant con- 
tribution to or at least a comparable one to that due to the "high energy" term 
oc |i?i3p|t/T-3p sin 2(^13. Accordingly, it is useful to know the values I-R12I and I-R13I which 
maximise the function 

FMRMRi'.\) = - 'vf"^"' . (23) 

The maximum of Fi(|i?i2|, |i?i3|) is obtained for [9] \R12\l\Rn\ = {Ami/ Amiy/^ = 2-4, 
and at the maximum we have F™^^ = 0.5 {Am\/ Amiy/^ = 1.2. At iRul/lRis] = 
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(Am^/AmQ)^/^, the corresponding function in the "high energy" term in er, 

Fs{\Ru\, \Ri3\) = , (24) 



=^ \Ri2? + \Rn? 
Am^ y I I I I 

takes the value 0.5, which is smaller only by a factor of 2 than its maximal possible value. 
The latter, however, takes place at I-R12I = 0, for which er = 62 = 0. 

The wash-out mass parameters in the case of interest are given by: 



rrir = ^JAm'Q\Rl2\'\Ur2r+^JAmi\Rl3\'\Ur3r 

+ 2(Am|Ami)^/^|i?i2| |i?i3|Re (e^^'^^^-'^^^f/.^af/.s) 



m2 



(25) 



Am| |/?i2|' + V Ami l^isl' - • (26) 



Note that the "high energy" phase difference ((/?i2 — v^is) in the expressions for rrir and 
m2 adds up to the Majorana phase (difference) 032/2, where 0:32 = (031 — «2i)- 

We will analyse next the combined effects of the "high energy" and "low energy" 
CP-violating phases on the generation of the baryon asymmetry. 



3.1 CP Violation Due to Majorana Phase in t/pMNS and i?— Phases 

Consider first the possibility that the baryon asymmetry \Yb\ is generated by the 
combined effect of CP-violation due to the Majorana phases in the PMNS matrix U and 
the phases (fi2 and (fis of the orthogonal matrix R. The Dirac phase 6 will be assumed to 
take a CP-conserving value: 6 = irk, k = 0,1,2, .... The CP- violating asymmetries and 
62 and the wash-out mass parameters and m2 depend in the case under study on the 
Majorana phase difference 032. We will be interested in the interplay between the effects 
of the phases 032, <^i2 and (^13 on the predicted value of \Yb\. 

In the case under consideration the asymmetry e,- can be written in the form: 



3 Ml J Ami r ^ ,9 

2 { ^3 Ur3 ^ sin 2(^13 



1 



Ami) 



^32. , / Am|\ 2 . ^32. 



.(27) 



where (/?23 = <^i2 + <^i3, the functions Fi and F3 are defined respectively by eqs. (23) 
and (24) and we have used the fact that for 6 = nk, (exp(— m32/2)[/*2^r3) = — (C12S23 ± 
■Si2C23Si3)c23Ci3 = —\U*2Ur3\. The asymmetry €2 can be obtained from eq. (27) by replac- 
ing |?7r3p with (1 — |t/T3p) and by changing the minus sign in front of the term oc Fi to 
plus sign. The expression for the baryon asymmetry has the form 

Yb = Y^ (Ahe + ^mix) , (28) 
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where 



i.-^i£LV:M^3xlO-f^Ur-™T7|. (29) 



37^* 16 7ri;2 VlO^GeVy I 5 x IQ-^ eV 



Ahe = Fs sin 2(^13 [ \Ur3f riiOMmr) + (1 - It/.s!') ^(O.Tlms) ] , (30) 
is the "high energy" term and 



2xi 



^Mix= - Fi|f/;2f/.3| [^(0.66m.)-r/(0.71m2) 



X 



1 



/ ~ ~ tt32 N / Am?) \ 2 ^ Q/32 ^ 

sin((^i2 + ifis + — ) + -—^ sin((^i2 + (^13 - — . 



(31) 



is the "mixed" term. Note that for the best fit value of S23 = 0.5, we have |t/r3p 



23""13 



0.5 ^ (1 - |f/r3^^/^ and therefore effectively Ahe oc (77(0. 66m^) + //(O.Tlma)). 
For ipi2 = kTi/2, <^i3 = A;'7r/2, k,k' = 0,1,2,..., one has Ahe = and we recover the 
expression for from [9] when the only source of CP- violation are the Majorana phases 
in the PMNS matrix U. One can have successful leptogenesis in this case typically for 
Ml ^ 4x10^° GeV and | sin 032 /2| ^ 0.1 [9]. The phase q;32 is also present in the expression 
for the neutrinoless double beta {{PP)ou-) decay effective Majorana mass corresponding 
to the NH spectrum (see, e.g. [30, 31]). 

Few more comments are in order. It follows from eq. (28) that the r and (e + n) 
CP-violating asymmetries generated by the "high energy" term always add up, while the 
r and (e + fi) asymmetries due to the "mixed" term tend to compensate each other. 
The contribution of the "mixed" term to Yb has the additional "suppression" factor 
{AuIq/ Am\y^^ = 0.42 in comparison to that due to the "high energy" term. For sin(<^i2+ 
<^i3 + 0^32/2) = 0, the "mixed" term |v4mix|, as can be shown, is smaller at least by the 
factor (Am|/Ami)^/2ci2/\/2 ^ 0.11 than the "high energy" term |Ahe|- Finally, the 
sign of is determined by the sign of sin2(^i3, while the sign of Auix depends on the 
signs of sin((^i2 + (^13 + 032/2) and (77(0. 66m^) - 77(0. 71m2)). 

The "high energy" term Ahe oc F3sin2<^i3 will be suppressed and will give a sub- 
dominant contribution in \Yb\ if either the phase of Rf^ is to a good approximation 
CP-conserving so that sin2(^i3 = 0, or I-R13I/I-R12I is sufficiently small. For sin2(^i3 = 
sin(9?23 + <^23) = and I-R13I, I-R12I 7^ 0, however, we also have sin((y923 — ^^23) = 0, implying 
that i?i2 -^13 TeaA, while RuRn is real or purely imaginary. This case has been 
studied in detail in [9, 16]. If, on the other hand, I-R13I = 0, we will have e,- = 62 = 0, 
and, as a consequence, Yb = 0. Thus, in order to have successful leptogenesis in the 
case of interest, the ratio I-R13I/I-R12I should not be too small, i.e. should be larger than 
approximately 0.05. 
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For the wash-out mass parameter rrir we have: 

= ^Aml\R^2\^\Ur2\^ + -^Ami li^isHf/^gp 

- 2 {Ami Ami)i/^ |i?i2| |i?i3| |f/;2f/r3| cos ((^12 - ^13 + ^) . (32) 
Thus, for given I-R12I and I-R13I, m,- satisfies the following inequalities: 

> ^MlU.,^{^-{^y''{l^0, (33) 

< y^|^,,nt/.p(l,Q)"|^|M)\ ,34) 

It follows from eq. (26) that the minimum (maximum) value of rhj. corresponds to the 
maximum (minimum) value of rh2- 

We have seen that for fixed Mi and given values of the neutrino oscillations parameters, 
the asymmetry Yb and the relative contributions to Yb of the "high energy" and the 
"mixed" terms depend on I-R12I, I-R13I and the Majorana phase 032, or equivalently, on 
the three phases ipu, <^i3 and 032- One of the constraints (^12 and (^13 should satisfy, as we 
have already indicated, is sgn(sin 2(^12) = —sgn(sin 2(^13) (see eq. (15)). It follows from 
eqs. (28) - (31) and (32) that ys(<^i2, <^i3; "32) = -yB{-<f 12, -<f 13] 4:7t ~ 032)- Therefore 
in what follows we shall analyse only the case of sin2(^i2 < 0, sin2(^i3 > 0. The results 
corresponding to sin2(^i2 > 0, sin2(^i3 < can be obtained from those we will derive 
using the indicated property of Yb- In what concerns the values of I-R12I and I-R13I, there 
are several possibilities leading to quite different physics results: 

i) I-R13I < I-R12I with |i?i2| < 1; ii) I-R12I < I-R13I with |i?i3| < 1; iii) |i?i2| > 1 or |i?i3| > 1. 

Consider first the case of I-R13I < I-R12I < 1- Under these conditions one always 
has I-R12I > l/\/2- As we have already indicated earlier, the asymmetry \Yb\ will be 
strongly suppressed if I-R13I/I-R12I <^ 0.05, so we will limit our discussion to values of 
I-R13I/I-R12I ~ 0.05. The results we obtain depend on whether |i?i3| < 0.5 or |_Ri3| > 0.5. 

For |i?i3| < 0.5, one should have |i?i2| > vO/75 = 0.87 in order for sin2<^i3 7^ 0. Let 
us set I-R12I = 1 for concreteness. In this case we get cos2(^i2 = 1 — 0.5|i?i3|'' ^ 0.97, 
|sin2(^i2| = |i?i3p < 0.25, cos2(^i3 = 0.5|i?i3p < 0.125, |sin2(^i3| = 1 - |i?i3|V8 > 
1 - 7.8 X 10-3. Thus, < (-^12) ^ 0.12 and (^13 = 7r/4. This implies that for ^32/2 = 
7r/4 we would have sin((^i2 + (^13 + 032/2) = 1, while if 032/2 = 37r/4 the "mixed" 
term will be strongly suppressed. It follows from these simple observations that the 
predictions for \Yb\ will exhibit a strong dependence on 0^2- For 032/2 = ir/A we also 
have cos((^i2 — <^i3 + 032/2) = cos (^12 = 1, and for any given |_Ri3| < 0.5, fhr will take to 
a good approximation its minimal value. 

Let us analyse the behavior of Auix and A^e for 032/2 = 7r/4 as I-R13I decreases 
starting from the value I-R13I = 0.5 (Fig. 1). In this analysis we will set sin^i3 = for 
simplicity and will use the best fit values of the other neutrino oscillation parameters 
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given in eq. (2); we will comment on the possible effects of sin ^13 having a value close to 
the existing upper limit later. It is useful also to recall that the efficiency factor r){X) i) 
has an absolute maximum a.t X = 10"^ eV, and at the maximum 77 = 7 x 10~^, and that 
ii) both for increasing X > 10~^ eV and decreasing X < 10"^ eV, f]{X) is a monotonically 
decreasing function of X. 

At |i?i3| = 0.5 we have: = 5.7 x 10"^ (weak wash-out), m2 = 2.1 x 10"^ > 
(strong wash-out), r/(0.66m^) = 4.2 x lO'^, and r/(0.71m2) = 6.8 x 10"^ < r/(0.66m^). The 
"mixed" term and the "high energy" term have opposite signs and Amix — —7.3 x 10~^ 
and Ahe — 1-40 x 10"^. Therefore the contribution of the "mixed" term in Yb has 
the effect of partially compensating the contribution of the "high energy" term, so that 
the sum (Amix + ^he) is approximately by a factor of 2 smaller than Ahe- As I-R13I 
decreases starting from 0.5, fhr, m2 and 77(0. 66mr) also decrease starting from the values 
given above. However, 77(0. 71m2) increases and at I-R13I = 0.41 one has 77(0.66772,-) = 
7/(0.71777,2). As a consequence, at I-R13I = 0.41, |Amix| goes through a deep minimum (at 
which it can even be 0) and is strongly suppressed. The "high energy" term |Ahe| just 
decreases somewhat as I-R13I changes from 0.50 to 0.41. At I-R13I = 0.41, the mixed term 
Amix changes sign: in the interval |i?i3| = (0.3 — 0.4) one has 7/(0.667/1,^) < 7/(0. 71m2) 
and, consequently Amix > 0. Thus, Ahe and Amix have the same sign and add up 
constructively in Yb- As I-R13I decreases below 0.41, rhr, m2 and 77(0. 66m^) continue to 
decrease, while 7/(0.7l7n2) continues to increase; Amix also increases rapidly, while Ahe 
decreases but rather slowly (Fig. 1). At I-R13I = (A777,Q/Am^)^/^ci2 = 0.35 we get fhr = 
and Amix has a local maximum. At this point we have Amix — Ahe = 2 x 10"^. As |i?i3| 
decreases further, 771^ and 7/(0.66m,r) increase, 7712 decreases, but 7/(0. 71m2) increases. 
As a consequence, Ahe also increases, while Amix diminishes. At I-R13I = 0.27 we have 
7/(0. 66m^) = 7/(0.71777,2) and Amix exhibits a second deep minimum, Amix — 0. At 
values of |i?i3| < 0.27 the inequality 7/(0. 66m,^) > 77(0.71772,2) holds and Amix is negative, 
Amix < 0. Therefore Ahe and Amix have opposite signs and their contributions to Yb tend 
to compensate each other. For decreasing |_Ri3| < 0.27, 7/(0. 66m^) and ^1(7/(0.66777,^) — 
7/(0.717772)) grow faster than 7/(0.717772) and ^3(7/(0.66777,^) -1- 7/(0.717772)), respectively. At 
iRi^l = 0.18, Ahe has a local maximum. However, one also has |Amix| — Ahe- As a 
consequence, Amix + Ahe — 0, i.e. the "high energy" and the "mixed" terms cancel 
each other and \Yb\ is strongly suppressed. This important feature of \Yb\ persists for 
values of 032/2 up to ~ n/2. The precise position of the considered deep minimum of 
\Yb\ depends on the value of sin^ 6^23 = S23 and, to less extent, on whether (5 = or vr if 
sin ^13 = Si3 has a value close to the existing upper limit. As an illustration, in Fig. 2 we 
show IY^AheI, I^b^mixI and 11^1 as functions of I-R13I for S23 = 0.64, S13 = 0.2 and 6 = 0. 
As is seen in the figure, for s^.^ = 0.64 we have Amix + Ahe — 0, and correspondingly 
Yb = 0, at |i?i3| = 0.30. Note that both |Y^Ahe| and |Y^Amix| have relatively large 
values at |i?i3| = 0.30 and thus each of the two terms separately could account for the 
observed value of Yb (Fig. 2). Nevertheless, the generated baryon asymmetry is strongly 
suppressed, \Yb\ = |^b(^he + Amix)| ^ 8.6 x 10"^^ and it is impossible to reproduce the 

^"This value is obtained as a solution of the equation 0.66mr/(8.25 x 10"'^ eV) = (0.71m2/(2 x 

10-^ ey))-i-i6. 
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measured value of Yb for Mi ^ 10^^ GeV. 

For I-R13I < 0.17, the "mixed" term is larger (in absolute value) than the "high energy" 
term, |Amix| > ^he; at |i?i3| = 0.10, for instance, we have |v4mix| — SAhe- Since the 
two terms have opposite signs, sgn(^Mix) = — sgn(^HE)5 the contributions of the "high 
energy" term in Yb compensates partially the contribution of the "mixed" term. 

We have studied also the dependence of the baryon asymmetry \Yb\ on the Majorana 
phase for given fixed values of I-R13I < I-R12I = 1 from the interval 0.1 ^ I-R13I ^ 0.5. 
This was done for three values of s^^ = 0.36; 0.50; 0.64 and two values of S13 = 0; 0.2. In 
the case of S13 = 0.2, the two CP-conserving values of the Dirac phase 6 were considered: 
6 = 0; 71. These results are illustrated in Figs. 3-5. As these figures indicate, the 
behavior of |Yb| as a function of 032 exhibits particularly interesting features when 
changes in the interval < 032 ^ vr. We note only that for S13 = 0.2 and given S23, we 
can get very different dependence of \Yb\ on 032 for the two values of 5 = 0; vr (Fig. 3), 
and that the dependence under discussion for, e.g. S23 = 0.50 can differ drastically from 
those for S23 = 0.36 and for S23 = 0.64 (Figs. 4-5). 

We can analyse in a similar manner the behavior of Auix, ^he and \Yb\ as functions 
|i?i3| in the interval 0.5 < I-R13I < 1.0. As in the preceding discussion we assume that 
I -Rial < I-R12I < 1.0 and fix again 032/2 = 7r/4. As can be easily verified, when I-R13I 
increases from 0.5 to 1.0 under the indicated conditions, i) Fi sin(<^i2 + (^13 + 032/2) 
changes from 1.14 to 0.60, ii) Fssin2(pi3 increases from 0.59 to 0.74, iii) rhr increases 
monotonically by a factor of ~ 20 from 5.7 x lO"'^ eV to 1.1 x 10~^ eV, and iv) m2 
increases only by a factor of ~ 2.3 from 2.1 x 10~^ eV to 4.8 x 10~^ eV. Correspondingly, 
the efficiency factor ?7(0.66mT-) first increases starting from the value 4.2 x 10"^, reaches 
a maximum 77(0. 66m^) = 6.8 x 10~^ at |i?i3| = 0.6 when 0.66777.^ = 1.1 x 10"^ eV, and 
then decreases monotonically to 1.52 x 10"^. In contrast, when |i?i3| changes from 0.5 
to 1.0, the efficiency factor 77(0. 71m2) only decreases monotonically by a factor of ~ 2.6 
approximately from 6.7 x 10~^ to 2.6 x 10~^. Thus, the asymmetry in the (e -|- /i) lepton 
charge is generated in the regime of strong wash-out, while the wash-out effects in the 
production of the asymmetry in r lepton charge change from weak to strong passing 
through a minimum. Clearly, the change of Auix and Ahe with |i?i3| is determined 
essentially by the behavior of 7/(0.66772,-). We also have ri{0.Q6rhr) > 7/(0.717772) in the 
case under discussion, implying that sgn(AMix) = — sgn(y4HE)- For the considered range 
of |i?i3| one typically has |Amix| — (0.5 — 0.6)Ahe, so there is a partial cancellation 
between the two terms Amix and Ahe in Yb (Fig. 1). 

It should be clear that Amix, ^he and \Yb\ will exhibit a different dependence on 
|i?i3| varying in the range 0.05 ^ I-Risl < I-R12I < 1 if ^32/2 differs significantly from 
7r/4. If 032/2 = 37r/4, for instance, we get typically |Amix| ^ I^heI- For |_Ri3| ^ 0.5 this 
essentially is due to the fact that sin((^i2 + <^i3 -|- 032/2) -C 1, while for 0.5 < I-R13I < 1, 
I-R12I = 1, it is a consequence of the fact that 7/(0.66777,-) and 7/(0.71?772) have rather close 
values: when I-R13I changes from 0.5 to 1.0, 7/(0.66777,-)— 7/(0.717772) decreases approximately 
from 7.6 x 10"'^ to 2.7 x 10~^; at the same time the sum 7/(0. 66m,-) + r/(0.71m2) changes 
from 3 X 10^^ to 10~^, remaining by a factor ~ 4 bigger than (r/(0.66?77,-) — 7/(0. 71m2)). 
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One can perform a similar analysis in the case of I-R12I > 1 or I-R13I > 1. The results 
we have obtained for I-R12I > 1 are illustrated in Fig. 6, which shows the dependence 
of IY^AheI, I^b^mixI and of IYbI on |_Ri3| for |_Ri2| = 1.2, 032/2 = 7r/4 and 5^3 = 0.5, 
Si3 = 0, Ml = 10^^ GeV. The figure exhibits some typical features, namely, the relevance 
of the "mixed" term in the region close to the minimal allowed value of I-R13I, i.e. for 
I-R13I ^ 1. If I-R12I > 1 (say |i?i2| = 1.2 as in Fig. 6), I-R13P can take values in the interval 
(I-R12P — 1) ^ l-^i3p ^ (I-R12P + 1)- When I-R13P changes from its minimal value to its 
maximal value, the phase 2<^i3, as can be easily shown, decreases from vr to 0, whereas 
2(fi2 changes from to (— tt), so that one always has sin2(^i2 < 0. Obviously, at I-R13P = 
(I-R12P — 1) we have Ahe = since sin2(^i3 = 0, while for 032/2 7^ nk, k = 0, 1,2, one 
finds, in general, Auix 7^ 0. For the value of 032/2 = 7r/4 (Fig. 6), for instance, we get 
^Mix — —3.9 X 10~^. The salient features of the behavior of Ahe and Amix as functions 
of \Ri3\ shown in Fig. 6, can be understood qualitatively from the behavior essentially 
of F3 sin(2(^i3)?7(0.66mT-) and of Fi sin((^i2 + (frs + 032/2)77(0.667711-): both quantities first 
grow relatively fast monotonically as I-R13I increases starting from its minimal value, but 
the former grows faster that the latter. Since i) Ahe increases faster than | Amix I starting 
from 0, ii) except at the extreme values of |i?i3| where Ahe = 0, we have Ahe > 0, 
iii) |v4mix| increases starting from a finite value but Amix < 0, there is always a value 
of I-R13I relatively close to its minimal value at which Ahe = |Amix|- Obviously, at 
this point the baryon asymmetry is strongly suppressed: Yb = Y^^Aue + Amix) = 
(Fig. 6). The behavior of Ahe and |Amix| as I-R13I increases beyond the point at which 
Yg = 0, is basically determined by 77(0. 66mr), which goes through a maximum and 
after that decreases monotonically. Let us note also that at certain value of I-R13I > 1, 
sin((^i2 + (fis + 032/2) can go through zero and changes sign. As a consequence, Amix 
also can change sign. 

As the results described above show, in the case of NH light neutrino mass spectrum 
and CP violation due the "low energy" Majorana phases in t/pMNS and "high energy" R- 
phases, the predicted baryon asymmetry can exhibit strong dependence on the Majorana 
phase 032 if the latter has a value in the interval < 032 < vr (sin2(^i2 < 0, sin2(^i3 > 0), 
or Sir < 032 < 4:71 (sin2(^i2 > 0, sin2<^i3 < 0). In the most extreme cases we can have 
either = or Yb compatible with the observations in a certain point of the relevant 
parameter space, depending on the value of O32. 

3.2 CP Violation Due to Dirac Phase in ?7pmns and i?— Phases 

Consider next the possibility that the CP-violation in "flavoured" leptogenesis is due to 
the Dirac phase 6 in the PMNS matrix U and to the "high energy" phases (^12 and (^13 of 
the matrix R. The Majorana phase O32 will be assumed to take a CP-conserving value: 
O32 = irk, k = 0, 1,2, .... The expression for the baryon asymmetry Yb also in this case 
can be cast in the form (28). The "high energy" term Ahe is the same as in the Majorana 
and R- matrix CP- violation case and is given by eq. (30). The "mixed" term has the 
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following form for arbitrary 032: 
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-^1 C23 ci3 [77(0. 66m^) - 77(0.717712)] 
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0^32 
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(36) 



The wash out mass parameter mr is given by 



Ami 1^121' \Ur2\^ + V Ami l^isl' Wrsl' - 2 (Am| Ami)'/' I^i2| l^^isl C23C13 



C12S23 cos (pi2 - v^is 



^32 



S12C23-S13 COS ipi2 - V^IS + 



"32 



(37) 



For e.g. ^32 = 2'iTk, k = 0,1,2,..., and <^i2,<^i3 = 0, ±7r, R12 and -R13 are real, 
^HE = 0, while in the "mixed" term only the part proportional to $mix non-zero, 
^Mix ^ ^Mix 0- The CP- violation in leptogenesis in this case is entirely due to the 
Dirac phase 6 in the PMNS matrix and we recover the results obtained in [9] . In particular, 
one can have successful leptogenesis for Mi ^ 5 x 10^^ GeV provided |si3sin(5| ^ 0.1 
For a32 = and R12R13 > (-R12-R13 < 0), the baryon asymmetry \Yb\ has a maximum 
at = 0.75, Rj^ ^ 0.25 {RI2 = 0.85, Rj^ ^ 0.15). 

We are interested in the regions of parameter space where there is a noticeable interplay 
between the CP violation due to the "high energy" phases (^12 and (^13 and the CP 
violation due to the Dirac phase 6. Since the CP-violation effects due to the Dirac phase 
are always suppressed by the relatively small experimentally allowed value of S13, the 
regions of interest would correspond to (^13 ~ 0, ±7r/2, where Ahe is also suppressed. 
The case of (^13 ~ 0, ±7r/2, corresponds to I-R13I taking values close to the boundaries: 
I-R13P ~ I I-R12P T 1 I- 

Note that the "mixed" term A^ix contains a piece which does not depend on the Dirac 
phase 6. This 5-independent piece is multiplied by C12S23 which is approximately at least 
by a factor 7 larger than the corresponding mixing angle factor S12C23S13 in the 5-dependent 
term $mix- ^^e region I-R13P ~ | I-R12P =F 1 I, we also have sin(<^i2 -|- <^i3 -|- a32/2) = 
for 0:32 = irk, and the (5-independent term in A^j^ will ^Iso be suppressed. 

^^Valucs of S13 <; 0.1 are within the range to be probed by "near" future experiments with reactor D^, 
[46, 47]. Future long baseline experiments will aim at measuring values of sin^ 6*13 as small as 10~"'-10~'^ 
and at constraining (or determining) 6 (see, e.g. [28]). 
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We have performed a detailed numerical analysis of this region of parameter space 
for CP- violating values of the Dirac phase 6 and a CP-conserving Majorana phase 032. 
The results of this analysis show that successful leptogenesis can still be realized for 
l-RisP ~ I — 1 I and I-R12I ~ Moreover, in the cases we have considered, the 

effects of the CP-violating Dirac phase are relevant in order to reproduce the observed 
value of the baryon asymmetry. 

In Fig. 7 we show \Yb\ as a function of 6 for |i?i2| = 1, S13 = 0.2, i) 032 = (left panel) 
and ii) 032 = vr (right panel). We choose |i?i3| close to its lower bound. In both the shown 
cases, there is a significant interference between the "high energy" and the "mixed" terms 
that can suppress or enhance the baryon asymmetry. The latter is controlled by the Dirac 
phase S. 

Thus, it follows from our analysis that if the Majorana phase 032 possesses a CP 
conserving value, there still will be regions in the parameter space where the effects 
of the CP-violating Dirac phase in the PMNS matrix can be significant in "flavoured" 
leptogenesis, even if CP-violation is due also to the "high energy" i?-matrix phases. 



4 Inverted Hierarchical Light Neutrino Mass Spec- 
trum 

We get very different results in the case of IH light neutrino mass spectrum. Now we have 
777.3 -C mi^2 — v^l Am^l = 0.05 eV. As in the analysis of NH mass spectrum, in what 
follows we neglect the effects of the lightest neutrino mass, 777,3. We also set I-R13I = 
{N3 decoupling [8]) for simplicity. The orthogonality condition for the elements of the R 
matrix of interest reads: Rfi + = 1. It leads to constraints similar to those given in 
eqs. (14) and (15). They permit to express the high energy CP-violating phases relevant 
for the discussion in this Section, ipu and <^i2, in terms of the absolute values of |-Rii| and 
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cos2v3i2 = 1^ — sin 2(^12 = =fV1 - cos^ 2(^i2 . (39) 
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The coefficients |-Rii| and I-R12I are constrained by conditions similar to those given in 
eqs. (18) and (19); they can formally be obtained from the latter by replacing I-R13I with 
We shall have also |/2iipsin2(^ii + |_Ri2p sin2(^i2 = 0. Using this condition one 
can write the CP violating asymmetry in the form: 
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For If II = /c7r/2, if 12 = k'7r/2, k,k' = 0,1,2,..., Ru and R12 are either real or purely 
imaginary and the expression for reduces to the one derived in [9] . As was noticed in [9] , 
the asymmetry will be strongly suppressed (rendering successful leptogenesis impossible) 
if sin((^ii + If 12) = 0, i.e. if RuRu is real; if, however, R11R12 is purely imaginary, 
then I sin((^ii + <fi2)\ = 1 and one can have successful leptogenesis with CP-violation due 
exclusively to the Majorana or Dirac phases in the PMNS matrix. 

It is easy to convince oneself using the expression for e,- given above and eq. (13) 
that in spite of the presence of "high energy" CP violation, the following relation holds 
in the case under discussion: €2 = — Cr- Thus, the baryon asymmetry Yb can be written 
as a function of e^- only, like in the case of the matrix R satisfying the CP-invariance 
constraints, i.e. having real and/or purely imaginary elements [9]: 

= Vi^f^r] ~ ^(^^2]] (41) 



37^, V V589 7 ' V589 
For the wash-out mass parameters we obtain: 



lAmil [\Rn\'\Url\' + \Rl2\'\Ur2\ 

+ 2 \Rn\ \Ri2\ Re {e'^^''-^''^ U;^Ur2) ] 



(42) 



m2 = \ \Aml\{\Ru\' + \R12Y) - (43) 



4.1 CP Violation Due to Majorana Phase in f/pMNS and i?— Phases 

We assume that the Dirac phase 5 = kir, k = 0, 1, 2, ... The baryon asymmetry takes the 
form: 

Yb = Y^{ A'i^ + )=Y^{CnE + Cuix ) ( V (0.66m.) - t] (0.71m2)) , (44) 

where Yg is given in eq. (29), A]^ oc Che is the "high energy" term, 0*^ Cmix is the 

"mixed" term, 

Che = Cn sin2(^n [\Uri\^ - \Ur2\^] , (45) 



Cmix — — Ci2|f/*if/r2| 



1 Ami (~ , ~ N • "21 „ . . - ^ - s "21 
9 Ta — 27 cos(v3ii + fi2) sm — + 2 sin(v9ii + fi2) cos — 

Zi /\Tn, * ^ ^ 



(46) 



and Cn = |i?ii|V(|^ii|' + 1^12!'), C12 = \RiiRi2\/{\Rii\^ + \Ri2\^). In deriving these 
expressions we have used the fact that for S = nk and the values of the neutrino mix- 
ing angles allowed by the data we have {exp{—ia2i/2)U*^Ur2) = — [■5i2Ci2(s23 — c^^sf^,) ± 
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S23C23Ci3Si3(si2 ~ c^2)] = —\U*iUr2\- Note that the first term in the square brackets in 
eq. (46) is suppressed by the factor 0.5AmQ/| Am^| = 0.016, which renders it practically 
negligible in what concerns the discussion that follows. Note also that in contrast to the 
case of NH spectrum, both the "high energy" and the "mixed" terms in Yb are propor- 
tional to the difference of the two relevant efficiency factors (77(0.66771,-) — 7/(0.71m2)). The 
wash-out mass parameter now reads: 



\Aml\ [\Rn\'\Uri\' + \Ru\''\Ur2\'' 

- 2\RiiRi2\\U:,Ur2\ COS ((^11 - <^12 + ^)] • (47) 

We see that both the "mixed" term A^ix ^^e wash-out mass parameter (and thus 
the efficiency factor in Yb) depend on the Majorana phase 021. The latter determines the 
range of possible values of the effective Majorana mass in neutrinoless double beta decay 
in the case IH light neutrino mass spectrum [31]. 

Consider next the contributions of the "mixed" and the "high energy" terms to Yb- 
We begin by noting that 

If/rlT - |f/r2p = (S?2 - cD-^L - 4S12C12S23C23S13COS5 

^ -0.20 - 0.92 si3 COS 5, (48) 
where we have used = 0-30 and S23 = 0.5. Similarly, we have 

\U*lUr2\ = kl2Cl2S23 - Si3S23C23(c?2 - S?2) C0S5| 

^ 0.50 |0.46 - 0.40 si3 COS 5|. (49) 

Obviously, for S13 = one obtains (|f/r2p — If^riP) = \U*iUr2\- However, if S13 = 0.2 
and 5 = 71, we get ([[/^ip - It/^aP) = -0.016, \U*-^Ur2\ = 0.27 and the "high energy" 
term will be strongly suppressed, being typically by more than an order of magnitude 
smaller (in absolute value) than the "mixed" term. Actually, as we are going to show, 
in a large region of the corresponding parameter space we have |Cmix| > |Che|, and 
therefore I^mixI > I^heI- certain the one indicated above, one can even get 

|C'mix| ^ |C*he|- In these cases the "high energy" term plays essentially no role in the 
generation of the baryon asymmetry. 

To be concrete, consider the dependence of the "high energy" and the "mixed" terms 
on I-R12I in the case when I-R12I < l-Riil — 1- Let us set first S13 = and let us work 
with sin 2(^11 < and sin2(^i2 > 0. In the region of I-R12I ~ 0.5, as can be easily 
verified, we have: |sin2(^ii| = I-R12P and therefore |Che| fx Gnl sin 2<^ii| oc |-Ri2p- At 
the same time <^i2 = vr/4 and, correspondingly, for, e.g. a2i/2 = 7r/4 we get: |Cmix| 
2G12I sin((^ii -|- (^12) cos a2i/2| cx; |i?i2|. Obviously, in the case under discussion one finds 
|C'mix| > |C'he| and the "mixed" term dominates over the "high energy" term, I^mixI > 
This interesting possibility is illustrated in Fig. 8. We note that Cmix and Che 



^^Onc can use the following property of the baryon asymmetry, yB(<^ii, (^12; Q;2i) 
-Yb(— (fill, —ipi2] 47r — a2i), to obtain results for sin2(pii > and sin2(pi2 < 0. 
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have opposite signs. The details of the behavior of I^mixI? I^heI of IYbI as functions 
of I-R12I, seen in the figure, can be understood following the behavior of | sin((^ii + (fu)], 
|sin2<^ii|, m,-, m2 and correspondingly of {ri{0.66rhr) — rj {0.7 1111,2)) ■ For example, the 
dominating maximum of I^mixI 1-^12 1 — 0-4 corresponds essentially to a maximum of 
{r]{OMmr) - r/(0.71m2)) at = 1.2 x 10"^ eV. For \Ri2\ ^ 0.5 and «2i/2 = ir/A we 
have cos(<^ii — (fi2 + tt2i/2) = 1 and for each given |i?i2| the wash-out mass parameter 
rhr has its minimal value (see eq. (47)). As I-R12I increases starting from the value of 
0.5, I sin 2(^11 1 also increases while | sin((^ii + <^i2)| and (77(0. 66m^) — ?7(0.71m2)) decrease. 
At I-R12I — 0.65, |(?7(0.66mr) — ?7(0.71m2))| goes through a minimum associated with rhr 
having a very small value, rhr — 0. At this minimum both I^mixI I^heI ^^^^ similar 
in magnitude and have relatively small values (Fig. 8). Correspondingly, \Yb\ is strongly 
suppressed and it is impossible to have successful leptogenesis. Actually, at I-R12I = 0.7 
one finds = 0. At this value of I-R12I the asymmetry Yb changes sign. As I-R12I 
increases beyond 0.7, |sin2<^ii| and (77(0. 66mT-) — ?7(0.71m2)) continue to increase, while 
I sin((^ii + <fi2) \ continues to decrease. As a consequence, at |i?i2| > 0.7 the "high energy" 
term is larger than the "mixed" term (Fig. 8) and the latter partially compensates the 
contribution of the former in Yg. At |i?i2| = 0.85, (//(O.GGm,-) — 77(0. 71m2)) has a second 
maximum. Now \A^\ dominates over I^mixI \^b\- 

Obviously, the results we get depend critically on the Majorana phase 0:21. If 021 = 
3n/2, for instance, \Yb\ will be strongly suppressed due to the factor {rj^O.GGfhr) — 
77(0.717722)) (strong wash-out regime). 

The possibility of strong suppression of the "high energy" term for S13 = 0.2 and 
(5 = vr is illustrated for I-R12I < l-Riil = 1 in Fig. 9. As the figure shows, in this case 
we have I^mkI ^ I^heI ^^e magnitude of the baryon asymmetry \Yb\ is determined 
entirely by the "mixed" term. A more detailed investigation of this interesting possibility 
is presented in [48] Obviously, the same conclusion is valid also for > 1, say 

|i?ii| = 1.2, if si3 = 0.2 and S = tt. 

Finally, in Fig. 10 we have presented graphically the results of a similar analysis 
performed for |-Rii| = 1.2, S13 = 0, and a2i/2 = it/A. The "mixed" term I^mixI ^^is 
case can be larger than, or comparable to, the "high energy" term \A^\ in a rather large 
interval of values of I-R12I which includes its minimal value, I-R12P ~ l-^iiP ~ 1; and in the 
"vicinity" of its maximal value, I-R12P ~ l-^iiP + 1- We can have a complete cancellation 
between the contributions of the "mixed" and "high energy" terms and Yb = at certain 

^•^The main differences between the analyses of the case of IH spectrum performed in the present 
article and in [48] are the following. Here we have investigated the case of R13 = 0, compatible with 
the hypothesis of decoupling of the heaviest RH Majorana neutrino N3. In [48] we have analysed the 
more general case of non-zero R13 and real i?x3- In the present paper we are interested in the general 
interplay between the contributions to the baryon asymmetry due to the "high energy" CP violating 
and that due to the "low energy" CP violating phases in the neutrino mixing matrix, while in [48] we 
have concentrated only on the cases in which the contribution to the baryon asymmetry due to the 
"high energy" CP violating phases (the i?-phases) is subdominant or strongly suppressed. Finally, in 
[48] we give a very detailed description of the regions of the leptogenesis parameter space, in which the 
contribution to the baryon asymmetry due to the "high energy" CP violating phases is subdominant or 
strongly suppressed. This was done both in the cases R13 ~ and of non-zero real i?x3- 
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value of |i?i2| (Fig. 10), although at the |i?i2| in question each of these two terms is 
sufficiently large to account for the observed value of the baryon asymmetry if the other 
term were not present. 



4.2 CP Violation Due to Dirac Phase in ?7pmns and i?— Phases 

This case corresponds to 021 = kn, k = 0, 1,2, .... The expression for the baryon asym- 
metry has the form given by eq. (44). The "high energy" term is the same as in the 
case of CP violation due to the Majorana and R phases and is determined by eq. (45). 
The form of the "mixed" term, to be denoted as C^j^ (and A^ix °^ ^mix) depends on 
whether 021 = {2k + l)7r or 021 = 2q7i, k,q = 0, 1, 2, 

Cmx = (-l)''2G'i2 sin(<^n + <^i2)c23S23Si3sin5, 021 = (2/c + l)7r , (50) 

Cmix - (-1)"^"^^ 26*12 sin((^ii + <^i2) [C12S12S23 + S13 {sj^ - C12) C23S23COS(5] ,^21 = 2qn .(51) 

Using the experimentally determined values of s^^ and 5^35 ^^e upper limit sf^ < 0.05, 
it is easy to convince oneself that the term involving the Dirac phase in C^y^, eq. (51), 
gives always a subdominant contribution to C^^. As a consequence, Cj^fx exhibits very 
weak dependence on 6 even for values of S13 close to the existing upper limit. The 
dominant CP-violating term is due to the i?-matrix. Therefore we will discuss only the 
case a2i = {2k + l)7r in what follows. 

For arbitrary 021, the wash-out mass parameter m,- reads: 



rrir = \/|AmA 



ci2|i?i2| - si2|/2ii|e*(^"-^^^+— ) j S23 

(52) 

( e,„l R-,„l -I- R,-, |p«(Vll-'^12 + - 



+ Si3C23e (^Si2|i?12| + Ci2|i?ll|e* 
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The maximum of the "mixed" term C^\-^ with respect to 6 occurs obviously for 6 = 
{2k' + 1)71/2, k' = 0,1,2,.... A more detailed analysis shows that for, e.g. 021 = ti", 
the absolute maximum of the efficiency term (?7(0.66m^) — 77(0. 71m2)) is reached in the 
case of 5 = 37r/2. The interplay between the CP-violation due to the Dirac Phase in U 
and -R— phases for 021 = vr and 6 = Sir /2 is illustrated in Figs. 11 and 12. In Fig. 11 
(Fig. 12) we show the dependence of |Y^A|^|, I^b^mix^I l^sl I -^12 1 i^i the case 

of I -Rill = 1 (I -Rill = 1.2). As Figs. 11 and 12 illustrate, there are substantial ranges 
of values of the relevant parameters for which the "mixed" term is of the order of, or 
exceeds, the "high energy" term. Typically the contributions of these two terms to Yb 
have opposite signs and tend to compensate each other. In certain points of the parameter 
space the compensation can be complete and we can have Yb = 0, although both I^^^heI 
and I^B^Mix^l "^^^ have relatively large values (see Figs. 11 and 12). In Fig. 13 we show 
the dependence of the baryon asymmetry on the Dirac phase 6 for 021 = ti", sis = 0.2 and 
one set of values of |-Rii| and I-R12I, at which the "mixed" term dominates in |Yb| when 
6 = 37r/2. As the figure indicates, Yb can exhibit very strong dependence on the Dirac 
phase 6 when the Majorana phase a2i takes the CP-conserving value a2i = n {Stt). 
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5 Conclusions 



In the present article we have investigated the interplay in "flavoured" leptogenesis bet- 
ween the "low energy" CP- violation, originating from the PMNS neutrino mixing matrix 
U and the "high energy" CP-violation which can be present in the matrix of neutrino 
Yukawa couplings, A, and can manifest itself only in "high" energy scale phenomena. We 
worked within the simplest type I see-saw theory with three heavy Majorana neutrinos 
Nj, j = 1, 2, 3, having a hierarchical mass spectrum with masses Mi <C M2^3. In the basis 
employed by us, the matrix of neutrino Yukawa couplings A is the only source of CP- 
violation in the lepton sector. In this analysis we used the "orthogonal" parametrisation 
of A, involving a complex orthogonal matrix R, R^R = RR^ = 1: A oc R^/rriU\ 
where M and m are diagonal matrices formed by the masses Mj > and > of Nj 
and of the light Majorana neutrinos z/^, j, k = 1,2, 3, This parametrisation proved rather 
convenient in the analysis performed in [9], which showed that the CP- violation necessary 
for a successful leptogenesis could be provided exclusively by the Majorana and/or Dirac 
physical phases in the neutrino mixing matrix U. It permitted to investigate also the 
combined effect of the CP-violation due to the neutrino mixing matrix U and the CP- 
violation due to the matrix R in the generation of the baryon asymmetry in "flavoured" 
leptogenesis. Throughout this study we used the terms "low energy" and "high energy" 
for the CP-violation originating respectively from the matrices U and R. The matrix R, as 
is well-known, does not affect the "low" energy neutrino mixing phenomenology. The two 
matrices U and R are, in general, independent. The source of the requisite CP-violation 
in "flavoured" leptogenesis can, in principle, be the matrix R, the PMNS matrix U, or 
both R and U. If the matrix R satisfies the general CP-invariance constraints (having real 
or purely imaginary elements [9]), while the PMNS matrix U does not satisfy these con- 
straints, we consider the CP-violation as originating from the neutrino mixing matrix U, 
i.e. from the Dirac and/or Majorana phases in U. This case has been studied in detail in 
[9]. If, however, the Dirac and Majorana phases in U take CP conserving values, while the 
matrix R, and the Yukawa couplings A do not satisfy the constraints following from the 
requirement of CP-invariance, the CP-violation will manifest itself only in "high" energy 
phenomena (like, e.g. leptogenesis) and will be due to the matrix R. As is well-known, 
one can have successful leptogenesis in this case as well. When neither U nor R satisfy the 
CP-invariance conditions, both U and R will be sources of CP-violation effects at "high" 
energies . In the present work we were primarily interested in this last possibility. 

In the case of hierarchical heavy Majorana neutrinos Mi <^ M2 <^ M3, the 

generated baryon asymmetry depends (linearly) on the mass of the lightest RH Ma- 
jorana neutrino A'^i, Mi, and on the elements -Ri^- of the matrix R, j = 1,2, 3, present in 
the neutrino Yukawa couplings of A^^i. The baryon asymmetry was assumed to be pro- 
duced in the "two-flavour regime" of leptogenesis, which is realised for Mi ^ 10^^ GeV. 
We have considered two types of light neutrino mass spectrum allowed by the existing 
data (see, e.g. [20]), namely, the normal hierarchical (NH), mi <^ 777,2 < 1^3, and the 
inverted hierarchical (IH), ma ^ rrii < m2. The lightest neutrinos mass in both cases 
was assumed to be negligibly small. Accordingly, in the case of NH (IH) spectrum we 
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have set also Rn = {Ris = 0), which is compatible with the hypothesis of decoupling 
of the heaviest RH Majorana neutrino N3 [8, 43, 45]. Under these conditions the orthog- 
onality of the i?- matrix implies the following constraint on the elements Rij of interest: 
RI2 + Ri3 = 1 + -R12 = !)• The phases of R12 and R13 {Rn and R12), '^12 and (^13 
{(fii and (^12), are the "high energy" CP- violating leptogenesis parameters. Using the 
orthogonality condition, <^i2 and (fis {ipu and ipu) can be expressed in terms of I-R12I 
and I -Rial (|-Rii| and |-Ri2|)- Thus, for fixed Mi and given neutrino oscillation parameters 
|Am^|, AitIq, sin^6'23, sin^ 6*12 and sin^6'i3, the baryon asymmetry in the cases we have 
studied depends, in general, on two "high energy" CP-violating {R-) phases and on one 
Majorana and one Dirac "low energy" CP-violating phases of the neutrino mixing matrix 
U. Under the conditions considered the Majorana phase which enters into the expression 
for Yb corresponding to NH (IH) spectrum is 032 = asi — 021 («2i)- 

Analyzing the possibility of NH spectrum and CP violation due to the Majorana 
phase 032 and the i?-phases (fu and (^13 (subsection 3.1), we have found that there exists 
a relatively large region of the relevant parameter space in which the predicted value of 
the baryon asymmetry exhibits a strong dependence on the Majorana phase 032 provided 
the latter lies in the interval < 0:32 < vr (if sin2(^i2 < 0, sin2(^i3 > 0), or Sir < 0:32 < 47r 
(when sin 2(^12 > 0, sin2<^i3 < 0). The regions typically correspond to 0.05 ^ I-R13I ^ 0.5, 
I-R13I < I-R12I < 1, and to |i?i2| > 1, l-RisP ~ — 1- Depending on the value of Q!32, 

we can have, e.g. either \Yb\ ^ 8.6 x 10"^^ or Yb compatible with the observations in 
the indicated regions. The effects of the "low energy" CP violation due to 032 can be 
non-negligible in leptogenesis also for 0.5 < I-R13I < I-R12I < 1 (Figs. 1, 2 and 6). In 
the regions of the parameter space where the Majorana phase effects are significant, the 
contributions to Yg due to the "high energy" CP-violation and that involving the "low 
energy" CP- violating phase 032 typically have opposite signs and tend to compensate each 
other. This mutual compensation can be complete and we can have Yb = for certain 
values of the relevant parameters, in spite of the fact that each of the two contributions 
can be sufficiently large to account by itself for the observed value of Yb- We have found 
also that in the regions of values of the parameters for which there is a significant interplay 
between the "high energy" and the "low energy" CP-violation, the predicted value of IYbI 
can exhibit strong dependence i) on the atmospheric neutrino mixing parameter sin^ 623 
when the latter is varied in the range (0.36 - 0.64) allowed by the data, and, ii) on whether 
the Dirac phase 5 = or tt, if the value of sin^ 6*13 is sufficiently large (Figs. 3-5). 

If in the case of NH spectrum the Majorana phase 032 has a CP-conserving value, 
a32 = 27ik, k = 0,1, 2, the "low energy" CP-violation is due only to the Dirac phase 6 
in the PMNS matrix U (subsection 3.2). A region in the corresponding parameter space 
where we have noticeable CP-violation effects due to 6 and successful leptogenesis still 
exists, but is very limited (Fig. 7). This is essentially a consequence of the suppression by 
the factor sin6'i3 < 0.22 of the CP-violation effects associated with 6. For CP-violating 
values of the "high energy" phases (pu and (^13 such that sin((^i2 + (fis) is significantly 
different from zero, the contribution of the term in Yb involving the Dirac phase is actually 
further suppressed. 

We obtained very different results for IH neutrino mass spectrum (Section 4). In 
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this case there are large regions of values of the corresponding parameters, for which the 
contribution to due to the "low energy" CP- violating Majorana phase a2i, or Dirac 
phase S (for a2i = {2k + l)vr), is comparable in magnitude, or exceeds, the purely "high 
energy" contribution in Y^, originating from CP- violation generated by the i?- matrix 
(Figs. 8 - 13). Moreover, in certain significant subregions of the indicated regions, the 
contribution to due to the "high energy" CP-violation is subdominant. We have found 
also that for (— sin6'i3 cos(5) ^0.1, the "high energy" term in is strongly suppressed 
by the factor (|f/rip — |f/r2n- The "high energy" phases (fu and (pi2 in this case can have 
large CP-violating values. Nevertheless, if the indicated inequality is fulfilled, the purely 
"high energy" contribution to Yb due to the CP-violating i?-phases would play practically 
no role in the generation of baryon asymmetry compatible with the observations. One 
would have successful leptogenesis in this case only if the requisite CP- violation is provided 
by the Majorana and/or Dirac phases in the neutrino mixing matrix. 

The results obtained in this study show that the CP- violation due to the "low energy" 
Majorana and Dirac phases in the neutrino mixing matrix can play a significant role 
in the production of baryon asymmetry compatible with the observation in "flavoured" 
leptogenesis even in the presence of "high energy" CP-violation generated by additional 
physical phases in the matrix of neutrino Yukawa couplings, e.g. by CP-violating phases 
in the complex orthogonal matrix R appearing in the "orthogonal parametrisation" of 
neutrino Yukawa couplings. 

Note Added. In the article [49] which appeared approximately 3 months after our 
paper, negative results regarding the possible effects of "low energy" CP violation in 
"flavoured" leptogenesis with hierarchical heavy (RH) and light Majorana neutrinos, when 
the "high energy" CP-violation is also present, were reported. It should be noted that the 
analysis performed in [49] differs substantially from the analysis performed here. In [49] 
the leptogenesis is considered in the framework of the SUSY extension of the Standard 
Model, more specifically, in the minimal Supergravity (MSUGRA) scenario with real 
boundary conditions, in which the dynamics responsible for supersymmetry breaking are 
flavour blind and all the lepton flavour and CP violation is controlled by the neutrino 
Yukawa couplings. The leptogenesis parameter space is constrained, in particular, by 
requiring that the /i — e-l-7 decay rate branching ratio, predicted in this scenario, satisfies 
BR{fi — »■ e -|- 7) > 10~^^. We work in the simpler non-SUSY version of leptogenesis. 
The difference between our results and those found in [49], as the authors of [49] also 
notice, may reflect the difference in the priors on the scanned leptogenesis parameters, for 
instance, the range in which the lightest RH neutrino mass Mi is varied. In any case, the 
aim of the analysis performed in [49] is different from the motivation and the resulting 
conclusions of our paper, as is clearly explained in the introduction of [49]. 
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Figure 1: The dependence of the "high energy" term IY^AheI (blue hne), the "mixed" 
term ly^^Mixl (green hne) and of the total baryon asymmetry \Yb\ (red line) on I-R13I in 
the case of NH spectrum, CP-violation due to the Majorana phases in U and i?-phases, 
= 7r/2, = 0.5, S13 = 0, |i?i2| = 1 and Mi = 10" GeV. 




IR13I 

Figure 2: The same as in Fig. 1, but for S23 = 0.64, S13 = 0.2 and 6 = 0. 
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Figure 3: The dependence of \Yb\ on the Majorana phase (difference) 032 in the case 
of NH spectrum, Majorana and R matrix CP-violation, = 0.5, Mi = 2 x 10^^ GeV, 

Rl2 = l, Rl3 = 0.19, i) Si3 

6 = 11 (blue line). 



(red line), ii) S13 = 0.2, 6 = (green line), iii) S13 = 0.2, 
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Figure 4: The dependence of IYbI on 032 in the case of NH spectrum, Majorana and R 
matrix CP-violation, \Ri2\ = 1, \Ri3\ = 0.51, Mi = 3.5 x 10^° GeV, sig = 0.5, S13 = 0.2 
and 5 = (tt) (red (green) line). 
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Figure 5: The same as in Fig. 4 but for = 0.64 
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Figure 6: The dependence of IY^AheI (blue hue), IY^AmixI (green hne) and of \Yb\ 
(red hne) on \Ri3\ in the case of NH spectrum, Majorana and R-matrix CP-violation, 
\Ri2\ = 1.2, 032/2 = 7r/4, = 0.5, S13 = and Mi = 10^^ GeV. 
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Figure 7: The dependence of \Yb\ on the Dirac phase 6 in the case of NH spectrum, Dirac 
and i?-matrix CP-violation, S13 = 0.2, = 1, Mi = 5 x 10" GeV and for i) 032 = 0, 
\Ri3\ = 0.16 (left panel) and ii) 032 = tt, iRnl = 0.12 (right panel). 




Figure 8: The dependence of the "high energy" term |Y^A[^| (blue line), the "mixed" 
term I^^^mixI (green line) and of the total baryon asymmetry \Yb\ (red line) on I-R12I in 
the case of IH spectrum, CP- violation due to the Majorana phase a2i in U and i?-phases, 
for 021 = 7r/2, = 1.0, S13 = and Mi = 10" GeV. 
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Figure 9: The same as in Fig. 8, but for S13 
I^b^heI (blue line) is strongly suppressed. 
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Figure 10: The dependence of the "high energy" term |1^J^he| (blue line), the "mixed" 
term I^^^mixI (green line) and of the total baryon asymmetry \Yb\ (red line) on I-R12I in 
the case of IH spectrum, CP-violation due to the Majorana phase 021 in U and i?-phases, 
for 021 = vr/2, \Ru\ = 1.2, S13 = and Mi = 10" GeV. 
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Figure 11: The dependence of the "high energy" term |5^^^he| (blue hne), the "mixed" 
term I^b^mixI (green hne) and of the total baryon asymmetry \Yb\ (red line) on I-R12I < 1 
in the case of IH spectrum, CP-violation due to the Dirac phase 6 in U and -R-phases, 
0.2, 5 = 37r/2, asi = rr, \Rn\ = 1-0 and Mi = lO^^ GeV. 
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Figure 12: The same as in Fig. 11 but for |-Rii| = 1.2. 
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Figure 13: The dependence of the baryon asymmetry \Yb\ on the Dirac phase 6 in the 
case of IH spectrum, Dirac and R phases CP-violation, S13 = 0.2, a2i = n, \Rii\ = 1, 
|i?i2| = 0.45 and Mi = 1.5 x 10^^ GeV. 
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